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A master equation was formulated to study intrinsic charge fluctuations of a grain in a plasma as 
ions and primary electrons are attached to the grain through collisional collection, and secondary 
electrons are emitted from the grain. Two different plasmas with Maxwellian and non-Maxwellian 
distributions were considered. The fluctuations could be bistable in either plasma when the sec¬ 
ondary electron emission is present, as two stable macrostates, associated with two stable roots of 
the charge net current, may exist. Metastablity of fluctuations, manifested by the passage of the 
grain charge between two macrostates, was shown to be possible. 


I. INTRODUCTION 

Various mechanisms, including ion and electron colli¬ 
sional collection and resulting secondary emission of ions 
or electrons, contribute to the charging of a dust grain in 
a plasma. Since the collision of plasma particles with the 
grain occurs at random times, the net electric charge pos¬ 
sessed by the grain fluctuates in time even if the plasma 
parameters such as temperature and number densities 
are fixed. This kind of fluctuations, which take place in 
systems with discrete particles, are known as intrinsic 
noise [T] . The intrinsic noise cannot be switched off as it 
is inherent in the actual physical mechanism, e.g., elec¬ 
tron or ion electron collision or emission in grain charg¬ 
ing mechanism, which is responsible for the evolution of 
the system. Intrinsic charge fluctuations refer to random 
variation of the grain charge by this intrinsic noise. 

Description of intrinsic charge fluctuations of grains 
was the subject of a number of studies PHIS]- Cui and 
Goree [2] studied the fluctuations through a Monte Carlo 
approach and concluded that they are more important for 
smaller grains. The grain charge is correlated with the 
grain size so this conclusion is consistent with that the 
net elementary charge Z possess ed by the grain should 
have fluctuations with Z^ms oc which suggested 

by Morflll et al. m- Cui and Goree [2] also showed 
that the fluctuating charge of small grains could expe¬ 
rience positive values. It is known that the grain mean 
charge at equilibrium is negative because in an average 
sense, the grain collects mobile electrons more than ions, 
as it approaches to an equilibrium charge. Matsoukas 
and Russell [3] proposed a one-step process master equa¬ 
tion [1] for the grain charge density function, then derived 
a Fokker-Planck equation for it and showed that if the 
condition e^/dTrcoR/c^Te ^ 1, where R is the radius of 
the grain and Tg is the electron temperature, is satisfied, 
the charge distribution at stationary states is Gaussian 
with an average and variance correlated to the ion and 
electron currents to the grain. Defining the system size 
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as U = AireoRkETe/e^, Shotorban [9] derived a Gaussian 
solution at non-stationary states for the Fokker-Planck 
equation formulated through the system size expansion 
of the master equation [1]. In the non-stationary state 
Gaussian solution, the rate of the mean grain charge 
correlates with the rate of charge to the net current. 
This mean equation is the macroscopic equation [T] of 
the grain charging system, and it is the same equation 
widely used for the grain charging with negligible fluc¬ 
tuations, which is the conservation of the grain charge. 
The rate of the grain charge variance is correlated with 
the currents and their derivatives evaluated at the charge 
mean. Shotorban m lately extended this model to in¬ 
clude multi-component plasmas where there are various 
kinds of singly- or multiply-charged negative or positive 
ions and showed that the grain charge distribution still 
follows Gaussianity when D is sufficiently large. In all 
the references discussed above, collisional collections of 
electrons and ions were the only mechanism of charging. 
Gordiets and Ferreira [6] obtained an analytical solution 
for the PDF at stationary states for a master equation 
that included the effect of the electron detachment, e.g., 
secondary electron emission (SEE), assuming that the 
grain charge does not experience positive values, i.e, a 
half-infinite range assumption Z = 0, —1, —2,— This 
kind stationary-state solution is unique for the master 
equation of a general one-step process with a half-infinite 
or finite range of the variable whereas it is not unique for 
a range consisting of all integers [T]. Later, Gordiets and 
Eerreira [8] formulated an improved version of the mas¬ 
ter equation that they had originally proposed m , re¬ 
laxed the half-infinite range assumption, and derived an 
approximate analytical solution for the PDE at the sta¬ 
tionary state. This approximation is not well justified for 
grains where the PDE varies substantially over a small 
range of charges. Khrapak et al. [7] studied the effects 
of thermionic emission and UV irradiation, separately, 
while electron collisional collection was present. They 
concluded that Zrms oc (Z) \ is valid for these situa¬ 
tions as well. Lately, Mishra and Misra m studied the 
fluctuations in multi-component plasma through a popu¬ 
lation balance equation resembling the master equation. 
They also included the influence of photoemission from 
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dust through irradiation by laser light in their study. It is 
noted that all above but Asgari et al [inun] used Markov 
approaches to describe grain charge fluctuations. 

The most well known effect of SEE on grain charging 
is a bifurcation phenomenon: two identical grains in an 
identical plasma environment may have two different sta¬ 
ble charge values, one positive and one negative [nnn]. 
Thus, a small variation in the parameters may cause a 
rapid change from one equilibrium charge to another. In¬ 
terestingly, Lai m showed that in the spacecraft charg¬ 
ing, the bifurcation phenomena caused by SEE could in¬ 
volve three stable equilibrium charge values. The exper¬ 
imental study of Walch et al m on charging of grains 
with energetic electrons that resulted in secondary elec¬ 
tron emission, showed the distribution of grain charge 
could be bimodal. However, they asserted that the lack of 
a unique value may be due to fluctuations in the plasma 
parameters or small differences in the grains. 

The current study is on the influence of SEE on grain 
charge intrinsic fluctuations with a focus on bistability. 
Bistability occurs in stochatic systems with two stable 
macrostates [Din]- The bistability of the grain charg¬ 
ing system is associated with the bifurication phenomena 
described above. The fluctuations in a bistable system 
may be metastable [T], where the fluctuations are at one 
macrostate for a while and at a random time, a passage to 
the other macrostate takes place and at a random time, 
the system returns to the first macrostate and this cy¬ 
cle continues. Whether grain charge fluctuations could 
be metastable is investigated in this work. In section [TT| 
first, a master equation describing the fluctuations of the 
grain charge in the presence of SEE mechanism is pre¬ 
sented, and then currents of ions, primary electrons and 
secondary electrons of a Maxwellian plasma and a non- 
Maxwellian plasma are shown. In section |ml results are 
shown and discussed. Conclusions are made in section 

M\ 


II. MATHEMATICAL FORMULATION 

Assuming that the charging of the grain undergoes a 
Markov process, the following master equation can be for¬ 
mulated for the probability density function of the grain 
charge P(Z, t): 

= (E - 1)MZ)I,{Z)P{Z) 

M-1 

+ ^ (E-" - 1) fn+l{Z)h{Z)P{Z) 

n=l 

+ (E-i - 1) Ii{Z)P{Z), (1) 

where E is an operator defined by W^g{Z) = g{Z ^ k) 
for any integer number /c, n indicates the number of sec¬ 
ondary electrons emitted from the grain upon the impact 
of one primary electron, M is the maximum number of 


secondary electrons that can be emitted, Ii{Z) and Ie{Z) 
are the currents of ions and primary electrons to the 
grain, respectively, and fn{Z) is the probability distri¬ 
bution of emission of n electrons in a single incident of 
a primary impact, i.e., the fraction of primary electrons 
that result in the emission of n secondary electrons in 
one single attachment incident. Hence, the rate of the 
attachment of the primary electrons that do not cause 
secondary emission is fo{Z)Ie{Z), and fn{Z)Ie{Z) indi¬ 
cate the rate of the attachment of the primary electrons 
that cause the emission of n secondary electrons in one 
incident. The jump process associated with eq. 0 is re¬ 
garded to that the attachment of a primary electron to 
the grain causes the emission of n secondary electrons; 
thus the net change of the grain charge is n — 1. In other 
words, Z{t)^ the charge of the grain at time t, jumps to 
Z(t) + n — 1. It is noted that the master equation of 
Gordiets and Eerreira [8] is a special case of eq. 0 with 
M = 3. Also, two following special cases of the master 
equation 0 regarded as one-step processes are worth 
noting: 

• M = 0, which corresponds to a case where no SEE 
occurs, i.e., fo{Z) = 1 and fn{Z) = 0 for n > 0. In 
this case, the second term on the right hand side 
of eq. 0 vanishes and the master equation of the 
grain charing is retrieved 13 E], 

• M = 1, which corresponds a case that at most one 
secondary electron is emitted so fo{Z) < 1 and 
h{Z) = l-fo{Z). 

Defining the system size U as a reference constant 
charge number and having changed the variable Z = 
where Z is modeled by a combination of 
a deterministic part (/)(t) scaled by U, and a random part 
^ scaled by a macroscopic equation associated with 
eq. 0 can be derived through the system size expansion 
method min]; 

f = <..(«, (2) 

where ai((/)) = J„(.) = /*{.) - /e(-) +4(-) is 

the net current to the grain, and Is{.) is the SEE current 
to the grain. A solution of the macroscopic eq. 0 is a 
time-dependent macrostate of the grain charging system 
while the solution of ai (0) = 0 is a stationary macrostate 
of the system [T] . 

Van Kampen [1] classifies the stable, bistable, and 
unstable stochastic systems through ai(0) in eq. 0 : 
A stochastic system is stable when a[{(l)) < 0 where 
a[{(l)) = dai/d(p; it is bistable when there are two sta¬ 
ble stationary macrostates, i.e., there are two solutions 
for ai(0) = 0 and at the vicinities of them, a[{(l)) < 0 
holds; and it is unstable when ai(0) > 0. The bistabil¬ 
ity classification here is in harmony with the bifurcation 
phenomenon based on the roots of the net current, which 
is reviewed in sec.|lj In other words, having neglected the 
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where for M > 0, p = n{Z)/M where n{Z) is given in 
eq. (§. Binomial distributions are used for the Monte 
Carlo modeling of SEE in the electron-surface collision 
[22] . Eor M = 1 in eq. (I^, fo{Z) = 1 — n{Z) and hence, 
fo{Z)Ie{Z)P{Z) = [/e(Z) - Is{Z)]P{Z). Eor M = 1, 
the summation terms in this equation are zero. It is 
noted that Ie{Z) — Is{Z) is the electron net current to 
the grain so the charging process here is modeled as each 
primary electron impact incident causing no or one sec¬ 
ondary electron emission at most. In eq. ([^, a sufficient 
condition for the positivity of fn{Z) is 1 —p > 0, which 
is equivalent to 


n{Z) < M. 


(6) 


FIG. 1. Dimensionless net current versus dimensionless 
grain charge z = Z/Q, where Q = ATveoRksTe/e^ and 
r = irR^ney/SksTe Imrie, for Ts/Te = 1.5 in a Maxwellian 
plasma m Case 1: Sm = 15 and EM/^kTe = 45.6; Case 
2: 6m = 14.85 and EM/^kT^ = 45.6; Case 3: 6m = 15 and 
EMl^kTe = 47. 


fluctuations of Z, i.e., ^ = 0, one obtains Z = Etcj) and 
eq. (§ is readily simplified to dZjdt = In{Z)^ which 
can be used to find the time evolution of Z. If an ini¬ 
tial Z is within the domain of attraction of a stable root 
of In{Z) = 0, then Z approaches it at the stationary 
state when there are no fluctuations. Each root is as¬ 
sociated with one stationary macrostate of the system 
and the system stability defined above can be likewise 
determined through the sign of I^iZ). When the fluctu¬ 
ations of the grain charge are taken into account, there 
is a probability for a fluctuation to carry the charge from 
the domain of attraction of one root to another. This 
situation is investigated in next section. 

The SEE current is correlated with Ie{Z) and fn{Z)^ 
and this correlation is found through the mean secondary 
electron yield defined by 


n(Z) = /,(Z)//e(Z). (3) 


This inequality sets the requirement for the minimum M 
according to variation of n{Z) versus Z. 

Eollowing Meyer-Vernet m, who investigated the 
bifurcation phenomena associated with SEE, both 
Maxwellian and non-Maxwellian plasmas are considered 
here: 

Eor Maxwellian plasmas, it can be shown HSlEa] 


Ie(z) = rx 


1 + ^ 


exp (g) 


Z>0, 
Z <0, 


( 7 ) 


Ii{Z) = Trii 




Z <0, 

(8) 

Z>0, 


where rii and Ti are the number density and temperature 
of ions, respectively. Also, Ti = TijT^^ fhi = m^/me, 

Ui = fli/Ue, 


In addition, n(Z) is correlated with fn{Z) through the 
definition of the mean n{Z) = Yl^=i'^fn{Z). Using 
these two equations and the normalization condition, i.e., 
fn{Z) = 1, one obtains 


Is{Z) = 


M 


l-foiZ) + J2in-l)fn{Z) 


n=2 


hiZ). (4) 


n = 


A-KeoRksTe 


( 9 ) 


r = -KR^ne 


SksTf, 


flujpeR 



( 10 ) 


Here, a binomial distribution is proposed for fn{Z): 


fn{Z) 


M 


P”(l-P) 


M—n 


( 5 ) 


where Xdc = electron Debye length 

and ujpe = nee‘^ f corrie is the electron plasma frequency. 

The SEE current is obtained by m 


n 
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FIG. 2. Probability density function of grain charge in a Maxwellian plasma; (a) Case 1 with R = 5nm; (b) Case 1 with 
R = 30nm; (c) Case 2 with R = 5nm; (d) Case 2 with R = 30nm; (e) Case 3 with R = 5nm; (C Case 3 with R = 30nm. 
Gaussian solutions are obtained by the system size expansion method lillS]. See the caption of fig.fpfor parameters associated 
with Case 1-3. 


ls{Z) 


3.76mT X 


^ exp (g) Fs 



Z > 0, 

Z <0, 


( 11 ) 


where 



Electrons in the non-Maxwellian plasma are assumed 
to have a bi-Maxwellian distribution. The electron cur¬ 
rent is obtained by adding an identical term, where Ue 
and Tg are to be replaced by uh and Th^ to the right 
side of eq. 0 and changing Ue to Ue — tih in the first 
term [17]. 


where B = ^AksTeZ/Q.E m‘ 
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FIG. 3. Grain charge variation in time in a Maxwellian plasma 
with SEE; (a) Gase 1 with R = 30nm and M = 3; (b) Gase 
2 with R = 5nm and M = 1; (c) Gase 3 with R = 5nm and 
M = 3. The dotted lines show the roots of the net current. 


III. RESULTS AND DISCUSSIONS 

Dimensionless net current is plotted against dimen¬ 
sionless grain potential in fig. for a Maxwellian 
plasma m- Seen in this figure is that when the SEE 
mechanism is lacking, the net current curve crosses the 
horizontal axis only at one point so there is only one root. 
The system is stable in this case as l!^{Z) < 0 for all val¬ 
ues of Z. This negativity is due to that the attaching 
electrons are more mobile than attaching ions. On the 
other hand, it is seen in the figure that when the SEE 
mechanism is present, the net current may have up to 
three roots, one negative and two positive (Cases 1 and 
2). The root at the middle is unstable whereas two oth¬ 
ers are stable so the grain charge fluctuations are bistable 
in Cases 1 and 2. In case 3, only one stable root exists 
and for all values of Z except the domain restricted be¬ 
tween local maxima and minima, the system is stable. It 
is noted that the SEE cases seen in fig. are different 
through small changes made in the SEE current parame¬ 
ters 6m or Em- A detailed study on the impact of these 
parameters and plasma parameters on the roots of net 
current can be found in Ref. nil 

The PDF of the grain charge obtained through a 
numerical solution of the master equation ^ for a 
Maxwellian plasma (Case 1-3 illustrated in fig.fff for two 
grain sizes R = 5 and R = 30nm are shown in fig. 

A prominent deviation from Gaussian distribution is ob¬ 
served for most of the cases. However, the PDF was 
verfied to be very close to Gaussian for M = 1 in figs. 



FIG. 4. Dimensionless net current versus dimensionless grain 
charge z = Z/D in a bi-Maxwellian plasma; Te — E — 25eV, 
6m — 3, EulksTe = 16, TsjTe = 1, Th iTe — 100, and 

M = 3 [II]. 


Hd -f) as compared to a Gaussian solution obtained by 
the system size expansion method I3II3]. For R = 5nm 
and M = 1, seen in fig. |^c), the PDF is bimodal, i.e., 
with two distinct local maxima, and for two other values 
of M, it is not. It is borne in mind that for all cases in 
figs[^a-d), there are two stable roots of the net current 
so they all are considered bistable according to the clas¬ 
sification in the previous section. The bimodal PDF is 
also seen in fig. [^b) for a larger grain with R = 30nm. 
However, when the same SEE parameters are used for 
smaller grain R = 5nm, no bimodal distribution is ob¬ 
served ^ee fig. If). No bimodal distribution is observed 
in fig. I^e-f) which is for the SEE cases with only one 
root of the net current. Although, in these two subfig¬ 
ures, the deviation of the distribution from Gaussianity 
is substantial for M = 2 and 3. 

Figure [^displays time histories of grain charges. The 
discrete stochastic method m, adapted from Gillespie’s 
algorithm [24l|25]), is utilized to simulate the grain charge 
fluctuations governed by the master equation Q. Time 
histories seen in figs. |^a,b), are for the bistable cases 
shown in fig. I^b) for M = 3 and fig. [^c) for M = 1, 
respectively. The fluctuations in these two cases are char¬ 
acterized by two distinct time scales: one associated with 
fluctuations around either of stable roots of the net cur¬ 
rent, i.e., charging macrostates, and the other associated 
with the spontaneous switches between them. A system 
with this behavior is called metastable [1]. A transition 
from the macrostate associated with the negative stable 
charge to the other macrostate is attributed to a sequence 
of incidents most of which increase the grain charge by 
one or two elementary charges. These incidents could be 
the attachment of an ion or the attachment of a primary 
electron that results in the emission of two or more of sec¬ 
ondary electrons. On the other hand, a transition from 
the macrostate associated with the positive stable root of 
the net current to the other macrostate is attributed to a 

































6 



FIG. 5. Probability density function of grain charge in a bi- 
Maxwellian plasma. See the caption of fig. I^for parameters. 




FIG. 6. Grain charge fluctuations in a bi-Maxwellian plasma 
with SEE, M = 3, and initial charge Z(0) = 0; (a) n^f/ne = 
0.035 and R = Inm; (b) nnlTie — 0.04 and R — Inm; (c) 
Tin!Tie — 0.04 and R — 3nm. 


sequence of incidents most of which are the attachments 
of a primary electron without emitting a secondary elec¬ 
tron. Fig. [^c) which corresponds to the PDF shown in 
fig- [^e) with M = 3 is not considered bistable as the net 
current in this case has only one root. 

Shown in fig. is the net current variation against the 
grain charge in a bi-Maxwellian plasma. For the shown 
SEE cases, there are two negative and one positive roots 
for the net current so the system is bistable in both SEE 
cases. The positive root is very close to the origin of the 
coordinates and the net current has a very sharp variation 
around this root. Although the curves of the SEE cases 
shown in this figure seem very similar, they are different 
as the negative roots in the case with = 0.035 are 

slightly closer to each other than the case with = 

0.04. The root of the net current in the No SEE case is 
somewhat far from the roots of the SEE cases. 

Figure displays the grain charge PDF in the stud¬ 


ied bi-Mawellian plasma for two grain sizes. All three 
cases shown in this figure are associated with the SEE 
cases in fig. which are bistable. A bimodal distribu¬ 
tion is observed for smaller grain with R = Inm at both 
= 0.035 and 0.04. Although, the difference be¬ 
tween these two values of is around %13, the bi¬ 

modal forms of their associated PDF’s are very different. 
The peak value of the PDF seen at around Z = 0 for 
= 0.035 is at least an order of magnitude larger 
than that for = 0.04. For this case, the value of 

the left peak is an order of magnitude larger that the 
right mode. For the grain with a larger radius R = 3nm, 
no bimodal behavior is observed. For this case, also, a 
Gaussian solution is obtained by the system size expan¬ 
sion with an initial condition (Z(0)) /O = —3. An ex¬ 
cellent agreement between the Gaussion solution and the 
master equation solution is observed. When (^(0)) = 0 
is used, the solution at the stationary state is a sharp 
Gaussian function at around Z = 0. Time history of the 
grain charge is shown in fig. with panels (a) and (b) as¬ 
sociated with solid- and dashed-line PDFs, respectively, 
in fig. An obvious metastability is observed for these 
two cases. 


IV. SUMMARY AND CONCLUSIONS 


A master equation was formulated to include the ef¬ 
fect of secondary electron emission in addition to colli- 
sional attachment of ions and electrons on the intrinsic 
charge fluctuations of a grain. Grain charging in both 
Maxwellian and non-Maxwellian plasmas were consid¬ 
ered. In both plasmas, the fluctuations could be bistable, 
as the system could have two stable macrostates. In the 
absence of SEE mechanism, the bistabillity is not possi¬ 
ble as the system always have a single macrostate. It was 
shown that if the system is bistable, the grain charge can 
be metastable. That is a situation where the fluctuations 
are characterized by two distinct time scales - one asso¬ 
ciated with the fluctuations at either macrostate and the 
other associated with the spontaneous transition between 
mactrostates. A switch from the macrostate associated 
with the negative stable root of the net current to the 
macrostate associate with the positive stable root of the 
net current is attributed to a sequence of incidents al¬ 
most all of which increase the grain charge by one or two 
elementary charge. On the other hand, a converse switch 
is attributed to a sequence of incidents most of which are 
the attachments of primary electrons without resulting 
in the emission of secondary electrons. 
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